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Abstract 

The local induction approximation (LIA) of the Biot-Savart law is often used for numerical and 
analytical investigations of vortex dynamics in the theory of superfluid turbulence. In this paper, 
using numerical simulation, some features of the LIA is considered. The temporal evolution of vortex 
loop spectrum is studied numerically. 
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1 Introduction. 

It is well known that chaotic vortex structures, appearing in volume of superfluid hehum under the 
particular conditions influence on hydrodynamic and thermodynamic properties of Hell. At this time, 
any explicit macroscopic equations to consider this influence are absent. In hydrodinamic approaches 
usually the Vinen equation is used. The system of hydrodynamic equations for superfluid turbulence, 
where the Vinen equation is incorporated, was obtained by some authors. (The details can be found 
for example in [1][2].) In turn, the Vinen equation is a phenomenological relation and it's validity is, 
in principle, unknown. From this point of view, for an adequate description of various phenomena in 
Hell, it is desirable to consider exact microscopic equations of vortex dynamics, based on the Biot- 
Savart law. But any analytical progress in this approach is rather difficult (practically zero at this time). 
Numerical investigations also have some problems. The Biot-Savart law implies the nonlocal interaction 
and therefore large computer resourses are demanded. Since the works [3] [4] to the present moment, the 
local induction approximation (LIA) of the Biot-Savart law is often used. But the validity of the LIA to 
describe vortex dynamics is the object of discussions until now. 

This report devoted to the numerical investigation of the LIA and present the continuation of the early 
analytical work [5] . In [5] was shown that the LIA can not describe a flux of amplitudes of harmonics 
across the spectrum in the one-dimesional Fourier-representation. This result seems rather unusual, 
taking into account that the LIA is a nonlinear differential equation. The aim of the paper is to check 
this with help of numerical modeling. 



2 The problem statement. 

The dynamical equation of the evolution of a vortex filament has the following LIA-form[3][4]: 



dt 



= 13s' X s" + vs' 
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Here s(^, t) is a radius vector of a vortex line point labeled by the variable ^, s' is the derivative on 
the parameter ^; t is time; the quantity /3 is the coefficient of the nonlinearity /3 = ^ log with the 
circulation k and the cutting parameters R (the external size, i.e. the averaged radius of curvature) and 
ro (the vortex core size). The coefficient of dissipation v appears in the eq.(l) when external counterflow 
is absent [3] [4] . Below one supposed that ^ is the arclength. 

Amplitudes of one-dimensional Fourier harmonics are defined as: 



The numerical calculations was executed for two cases. Firstly, the spectral evolution of some initial 
configuration was investigated. At that it was supposed i/ = in the eq.(l) and thus only nonlinear term 
was presented in the dynamical equation. The evolution of harmonic amplitudes was defined by this term 
only. Secondly, a random force term of the Gaussian type was added to the right hand side of eq.(l). 
Thus eq.(l) was modified as: 



3 The details of the numerical procedure 

The technical detail: the figures, mentioned below are available in the individual files. For example the 
file Fig.l.gif corresponds to Fig.l etc. 

In this chapter the full details of the numerical procedures are described. A reader, not interesting 
this theme, may omit this section. 

A vortex line was presented as a number of points in the coordinate space {x,y,z). For example, 
the ring with radius 1 contains 1200-5000 points, situated on the equal distance (As) from each other. 
The fourth order Runge-Kutta method was used for modeling the temporal evolution. The time step 
was defined as = (As)^/C. Where As is the minimal for the entire configuration space step, C is 
a constant. In the work [6] one maintain that this method is stable when C > \/2. The experimental 
testing in our case gives the result C ~ 500 — 1000 if the vortex line is strongly intricated. The criterion 
of the explicity was the conservation of the full length and the integral of the curvature square. (This 
values are some of the LIA invariants.) For example, this integral was equal 38222.2167961644 and it was 
covserved with the explicity ±0.0000001 during the program ran, the total length was conserved down to 
12 decimal signs. Besides to test the explicity, the iterative Crank-Nicolson type numerical scheme was 
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developed. This scheme was similar to the first one, described in the work [7]. But the distinction was 
the following: the equation for s was used, instead of the equation for a tangent vector I as in [7] . (There 
method is really unusable for the required explicity.) The criterion of the iteration convergence was 
about lO^^'^. Let's note that the explicity of the floating number representation in a standard computer 
is about 19 characters (the format "double"). For example, the same initial configurations (it is shown at 
the Fig.l on the (x,y) plane) were calculated with these methods till 192547 and 1174140 time steps, but 
for the approximately same physical time. One can see from Fig. 2 and Fig. 3 (the configurations on the 
(x, y)-plane, corresponding to these two calculations), the results are absolutely identical. The spectral 
distributions are also identical. 

To create a strongly intricated initial configuration, the following procedure was used. Usual ring, 
with radius 1 was created at first. After that, sinusoidal disturbances were added to this configuration. 
The disturbances were added step by step, and after each step the procedure of the rescaling of the length 
was applied, besides the coordinates of points were corrected to place them uniformly. In this way (step 
by step) some harmonics were excited. 

When the randomized model was studied, the additional /(i,^) was incorporated into the numerical 
method. At that this value was modeled in the following way: 



Here (-j^) is a one-dimensional wave vector kn with number n, Lt is a total length of a vortex line 
(recalculated at each time step), ^ is a current distance along the line, to clarify y see the rel.(4), a 
denotes a spatial component, is a random amplitude, defined by a random-number generator, (j)^ is a 
random phase. The last two values were generated at each time step. 

Besides, a random disturbance also was used to create an initial configuration. At that a bunch of 
harmonics was excited, and then the random force was switched off. 

During the was program running, one was possible to observe the changing of various parameters and 
the spectral picture on the monitor. To develop the program, Borland C+-|-Builder 5 tool was used. 

4 The results, discussion and conclusions 

Thus, in the previous paper [5] was shown that the LIA don't describe a flux of harmonic amplitudes 
across the spectrum. The aim of the numerical experiment was to confirm this result. The first part of 
the experiment looks as following. One was created an initial strongly intricated configuration with some 
spectral structure. After that, the vortex filament evolved accordingly to the eq.(l) (under the condition 
u = Q). The reconecting processes was omitted from the consideration, i.e. parts of a vortex line freely 
passed across each other when they intersected. The changing of the spectral picture was studed. The 
initial spectral distribution ploted on the Fig. 4 from n = 1 till n = 150 (there and below a red line 
correspond to an initial state, black line denotes a current distribution). A wave vector fc„ = 2-Kn/Lt, 
where Lt is the total length of the filament at the moment t. In coordinate spase this configuration 




n=nhigh 
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plotted on the Fig.l ((x,y)-plane). The successive steps of the spectral evolution are shown at the Fig.5- 
Fig.8. Finally, the spectral distribution looks as at the Fig. 9. At later time a stationary state was formed 
and any changes were only inside the area, marked by blue line. In the stationary state, the average 
size of the vortex tangle ceased to grow and only slight pulses occurcd. Inside of the marked area some 
peaks appeared, which dassapeared after a time. The general spectral distribution was not changed. The 
harmonic excitation in the low region (about the first harmonic) was absent just as the excitation in the 
high area. Thus, after the stationary state was formed, an spectrum changes was ceased. The initial 
spectral distribution was changed cardinally. One can seem that this point is not in agreement with the 
result of [5]. But the procedure, used in [5] implicity implies just stationary state. As well known, a 
standard situation look as following. A wave packet spreads to infinity when a nonlinear, nondissipative 
equation is used. 

The stability of the harmonic distribution, perhaps connected with the infinite number of the invari- 
ants, produced by the LIA. For example, total length L = const, J s" = const, etc. High harmonics 
can not excitate becaiise it means the growth of the ciirvature and, hence, of the curvature square, this is 
impossible. Increasing of low harmonics lead to the growth of a loop size, this is also impossible because 
-L is a constant and so the total curvature must decrease in this case. Thus some transfer across the 
spectrum is certanly absent if a vortex configuration achieved a stationary state. 

In the second part of the numerical experiment the temporary evolution of a detached narrow harmonic 
bunch was studied. The harmonics in the environment of this bunch was equal to zero. The initial 
configuration presented a ring. Some scores of time steps this ring imdergo the random distmbances of 
the kind: f"{t,^) = X]i="^''','' st*^'; " sin{ ^™^ + 0"). The amplitudes c" and the phases 0" possessed 
random values at each time steps. (Their values was defined by the random number generator.) As the 
result the initial spectral configuration was obtained (Fig. 10). (At the Fig. 10 the stripe of the spectrum 
is plotted from harmonic number n = 70 up to number n ~ 120.) Later the external force was switched 
off and the evolution under the LIA nonlinearity only was investigated. 

The result is following. There are no any additional harmonics exciting in the system, on condition 
that they equal to zero initially. This result was checked thoroughly. The boundaries of the excited area 
was thoroughly examined after about 800 thousands time steps. They were examined at the enlarge scale 
factor 10® (in the vertical direction) relative to the amlitude of the stripe. As it shown at the Fig. 11 where 
the boundaries of the excited harmonic bunch plotted (scale factor 10®), there are no any disturbances 
beyong the initial stripe. At the same time the spectral distribution inside of the perturbed area was 
varied. This change represented the periodical process (the period is about 20 thousands time steps). 
The steps of the changes is shown on the Fig.l2-Fig.l8. Obviously, the spectral distribution strictly 
return to it's initial state. The period of the return depends on the width of an excited strip and also 
on current numbers of the harmonics. Note again, in a standard nonlinear task we have a wave packet 
is smearing and tending to it's stationary spectral distribution. Note, however, after some millions time 
steps, distortion of the spectral stripe was observed. Apparently it is produced not strictly zero initial 
background. It was tested for the smooth ring. Obviously, in the later case only the first harmonic must 
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be present. But other harmonics had amphtudes of order 10~^ (relative to the amphtudes, plotted on 
the Fig. 12 - Fig. 18). Because, any nimcrical method has inadvertent errors, backgroud amplitudes was 
conserved within about 1 percent. At the Fig. 19 (scale factor 5 * 10^) shown the divergence from zero for 
the initial (red) and current (black) backgrounds. Besides, the graphic on the Fig. 11 is wider the on the 
Fig. 12 - Fig. 18 because the greater magnification in the vertical direction allows to see small amplitudes 
which are lost details. 

Also a case when the external force not switched off was considered. The force was during all time 
when the programm was running. At that, of course, i> ^ in the cq.(3). The result is the same. After 
the quasi-equilibrium between pumping and dissipation was formed (i.e. when the tendency for changing 
the total length and the total curvature vanished), any new harmonics not appeared in the system. 

Thus two conclusions for sure following for the LIA. 

1. There are no transfer across the spectrum when a vortex system is in the stationary state. 

2. There are only initially excited harmonics remain in a vortex system. New harmonics never excite. 
All redistributions of hamronic amplitudes take place between excited harmonics. 

And the last remark. It can be suppose that the question: " whether can LIA approximately describe 
the superflud turbulence" is rather make no sense until we don't know exactly the role of the processes 
of reconnections in dynamics of vortex loops. But this problem is too much difficult. 

The work was carried out under support of INTAS (grant N 2001-0618) and RFBR (grant N 03-02- 
16179). 
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